Ground-state structure of KNbOa/KTaOa superlattices: 
Array of nearly independent ferroelectrically ordered planes 
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Phonon spectra of the paraelectric phase, structure and distribution of polarization in the polar 
ground state, and energies of ferroelectrically and antiferroelectrically ordered phases were calculated 
for stress-free (KNb03)i(KTa03)n superlattices (n = 1-7) from first principles within the density- 
functional theory. Quasi-two-dimensional ferroelectric state with the polarization localized in the 
KNbOs layer was revealed to be the ground state for the superlattices with n > 2. The interaction 
between these polarized layers decreases exponentially with increasing the interlayer distance, and 
the ground-state structure transforms to an array of nearly independent ferroelectrically ordered 
planes at large n. 

PACS numbers: 68.65.Cd, 63.22.-m, 77.84.Ek, 81.05. Zx 
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I. INTRODUCTION 

Ferroelectric superlattices (SLs) have drawn much at- 
tention in recent years. It appeared that many physi- 
cal properties of these artificial structures, such as the 
Curie temperature, spontaneous polarization, dielectric 
constant, nonlinear dielectric and optical susceptibilities, 
exceed considerably the properties of bulk crystals and 
thin films of the corresponding solid solutions. 

Ferroelectric phenomena in KTai_a;Nba;03 solid solu- 
tions and KNbOa/KTaOs superlattices have been stud- 
ied for many years. The experiment shows that in solid 
solutions a polar phase appears at x > Xc = 0.008,^ but 
its properties in the vicinity of critical concentration Xc 
are modified by manifestation of quantum effects and by 
local disorder (the dipolar glass). To explain the pecu- 
liar properties of the solid solutions, a number of models 
has been proposed^i^ including a hypothesis about off- 
centering of Nb atoms in KTaOa^r— However, as was 
shown in Ref. the single-well shape of the local po- 
tential for Nb atom in the solid solution excludes its off- 
centering. 

Ferroelectric properties of KNbOa/KTaOa superlat- 
tices have been studied both experimentallyiSrJ^ and 
theoretically^^"— The ferroelectric properties of the SLs 
grown on KTaOs substrates were calculated assuming 
that the polar ground-state structure is tetragonal. The 
signs of antiferroelectricity in these SLs observed exper- 
imentallyi^ were confirmed by observation of unstable 
phonon at X point of the Brillouin zone in first-principles 
calculation^^ At the same time, our recent calculations 
of the ground state for both stress-free and substrate- 
supported superlattices with n = 1 and 3 (Ref. HO) re- 
vealed that the ground-state structure of KNbOa /KTaOa 
superlattices is monoclinic or orthorhombic, and is not 
tetragonal. Along with the ferroelectric instability, three 
unstable phonon modes that are not typical neither for 
KNbOs nor for KTaOs were observed at X, Z, and R 
points of the Brillouin zone in the paraelectric phase of 
these SLs. In contrast to other instabilities found in 



ten SLs studied in Ref. 20;, the origin of these instabil- 
ities was not clear. The calculations2i of the phonon 
frequencies at the F point for the paraelectric phase of 
(KNb03)i(KTa03)7 superlattice with the lattice param- 
eter equal to that of KTaOa revealed a weak ferroelectric 
instability of this phase with the polarization parallel to 
the layer plane. However, the distribution of polarization 
in the low-symmetry polar phase was not analyzed. 

The purpose of this paper was to study the de- 
velopment of the ferroelectric properties in stress-free 
(KNb03)i(KTa03)„ superlattices when increasing their 
period L = n + 1. We show that quasi- two-dimensional 
(2D) ferroelectric state with polarization localized in the 
KNbOa layer is the ground-state structure for the super- 
lattices with n > 2 and it transforms to an array of nearly 
independent ferroelectrically ordered planes at large n. 



II. CALCULATION DETAILS 

The (KNb03)i(KTa03)„ superlattices considered in 
this work are periodic structures grown in [001] direction 
and consisted of the KNbOa layer with a thickness of one 
unit cell and the KTaOa layer with a thickness of n unit 
cells {1 < n < 7). The modeling of these structures was 
performed using supercells with a size of lxlx(n+ 1) 
unit cells. 

The calculations were performed within the first- 
principles density-functional theory (DFT) with pseu- 
dopotentials and a plane-wave basis set as implemented 
in ABINIT software.— The local density approximation 
(LDA)— for the exchange-correlation functional was 
used. Optimized separable nonlocal pseudopotentials^i 
were constructed using the OPIUM software^S to improve 
the transferability of pseudopotentials, the local poten- 
tial correction's was added. Parameters used for con- 
struction of pseudopotentials and other details of calcu- 
lations are given in Ref. 2(X. The plane-wave cut-off en- 
ergy was 40 Ha (1088 eV). The integration over the Bril- 
louin zone was performed on a 8x8x4 Monkhorst-Pack 
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TABLE L Frequencies of unstable phonon modes at high- 
symmetry points of the Brillouin zone in the PA/mmm phase 
of (KNbOs)! (KTa03)n superlattices with n = 1,2, and 3. 



FIG. 1: Frequencies of lowest A2u and Eu polar optical modes 
in the paraelectric Pi/mmm phase as a function of the thick- 
ness of the KTaOs layer in (KNb03)i(KTa03)„ superlattices. 
Imaginary frequencies of unstable modes are plotted as nega- 
tive values. 



mesh for SLs with n = 1 and 2, and on a 8x8x2 one for 
SLs with n = 3-7. The relaxation of atomic positions 
and of the unit cell parameters was stopped when the 
Hellmann-Feynman forces were below 5 ■ 10~^ Ha/Bohr 
(0.25 meV/A). The phonon spectra were calculated using 
the density-functional perturbation theory. The sponta- 
neous polarization was calculated using the Berry phase 
methodi^l 

The structures of KNbOa and KTaOa calculated using 
this approach agree well with experiment. The calcu- 
lated lattice parameters for cubic KNbOa and KTaOa are 
3.983 and 3.937 A, they differ from the experimental val- 
ues (4.016 and 3.980 A) by 0.81 and 1.07%, respectively 
(small underestimation of the lattice parameters is typi- 
cal for LDA approximation). The calculated c/a ratio of 
1.0197 for tetragonal KNbOa is close to the experimental 
value (1.0165). The ground-state structure is rhombohe- 
dral {R3m) for KNbOa and cubic (PmSm) for KTaOs, in 
agreement with experiment. The calculated spontaneous 
polarization of 0.372 C/m^ for tetragonal KNbOsJs close 
to experimental values of 0.37-0.39 C/m^ (Ref. [28[ ). 



III. RESULTS 

Calculated frequencies of two lowest and £"„ 

ferroelectric modes at the F point of the Brillouin 
zone in the paraelectric PA/mmm phase of stress- free 
(KNb03)i(KTa03)„ superlattices with n = 1-7 are pre- 
sented in Fig. [T] In SL with n = 1 both the modes, 
which result from unstable TO phonon of the perovskite 
structure, are unstable. With increasing n, the 
mode becomes stable, and there remains one unstable 
Eu phonon mode, whose frequency is weakly dependent 
on the KTaOa layer thickness. 
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When searching for an equilibrium structure, all pos- 
sible structure distortions were tested. For SLs with 
n = 1, 2, and 3, in which two unstable modes appeared in 
the paraelectric phase, monoclinic phases with Cm and 
Pm space groups (polarized along [xxz] and [xOz] direc- 
tions) as well as tetragonal PAmrn and two orthorhom- 
bic Amm2 and Pmm2 phases were considered. For SLs 
with n > 3 only the Eu mode is unstable, and so only 
orthorhombic A'mm2 and Pmm2 phases polarized along 
[xxQ] and [xQO] directions were possible. The calcula- 
tions revealed that for SLs with n = 1 the ground-state 
structure is Cm^ whereas for other SLs — in spite of the 
existence of unstable A2u mode in the paraelectric phase 
of SLs with n = 2 and 3 — the ground-state structure is 
Amm2. 

Calculation of phonon frequencies at high-symmetry 
points of the Brillouin zone for paraelectric PA/mmm 
phase of SLs with n < 3 reveals, in addition to the fer- 
roelectric instability, unstable phonons at Z, X, and R 
points (Table HI). However, among all possible distorted 
structures the polar Cm and Amm2 phases had the low- 
est total energy, and all phonons at all points of the Bril- 
louin zone were stable in these phases. 

As follows from Table HIl the energy gain resulting from 
the ferroelectric ordering does not vanish with increasing 
n. This fact and the existence of unstable Eu phonon in 
the paraelectric phase give evidence for stability of the 
ferroelectric state in KNbOa layers of minimal thickness 
(one unit cell) located between layers of KTaOa of arbi- 
trary thickness. 

The polarization profiles for the superlattices were 
calculated using the approximate formula Psa — 



(e/f7)E.m.Z* 



where P^a is the a-component of 



polarization (a — x, y, z), VL is the unit cell volume of the 
layer under consideration, Uip is the /3-component of dis- 
placement of the ith atom from its position in the para- 
electric phase, Z* is the Born effective charge tensor of 
this atom,^^ and Wi is the weight equal to 1 for Nb(Ta) 
and O atoms lying in the Nb(Ta)-0 layer under consid- 
eration, 1/2 for K and O atoms lying in two nearest K-0 
planes, and for other atoms. The dependence of |Ps| on 
the layer number in SLs with ti = 3, 5, and 7 are given 
in Fig. [2l It is seen that in the region between KNbOa 
layers the polarization decays exponentially, with a char- 
acteristic length of ~3 A. This result indicates a strong 
localization of polarization in the potassium niobate layer 
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TABLE II: The total energy difference between paraelectric and ground-state structures, the lattice parameters, and the 
Berry-phase polarization for the ground-state structure of (KNb03)i(KTa03),i superlattices with different n. 



n 1 2 3 4 5 7__ 

AE (meV) 11.738 7.532 5.898 5.365 5.090 4.625 

a = 6 (A) 3.9656 3.9567 3.9508 3.9478 3.9459 3.9437 

c (A) 7.9235 11.8500 15.7908 19.7298 23.6680 31.5436 

7 (deg.) 89.9581 89.9746 89.9866 89.9913 89.9936 89.9958 

Q = /3 (deg.) 89.9701 90.0 90.0 90.0 90.0 90.0 

\Ps\ (C/m^) 0.2515 0.1192 0.0822 0.0612 0.0486 0.0338 




FIG. 2: The polarization profiles for (KNb03)i(KTa03)n su- 
perlattices with n = 3, 5, and 7. The number A'' = corre- 
sponds to the KNbOs layer. 



TABLE III: The energy difference between ferroelec- 
trically and antiferroelectrically ordered structures in 
[(KNb03)i(KTa03)„]2 superlattices (per Nb atom). 
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1 


2 


3 


2Wint (meV) 


4.282 


1.089 


0.278 



and proves the quiasi-two-dimensional character offerro- 
electricity in considered superlattices. The polarization 
in the KNbOa layer decreases monotonically from 0.279 
to 0.257 C/m^ with increasing n from 3 to 7 (for com- 
parison, for bulk orthorhombic KNbOa our calculations 
give \Ps\ = 0.418 C/m^). The total polarization of the 
supercells calculated using the above formula differs from 
the result of correct Berry-phase calculations (Table 
by only 3%. 

In order to determine the interaction energy Wint be- 
tween two polarized KNbOs layers separated by KTaOs 
layer, the total energies of [(KNb03)i(KTa03)„]2 super- 
lattices of doubled period with ferroelectric and antifer- 
roelectric ordering of polarization in neighboring KNbOs 
layers were calculated. The energy difference between 



the two ordered structures (per Nb atom) is given in Ta- 
ble mil As follows from these data, the energy differ- 
ence decreases exponentially with increasing thickness of 
KTaOa layer with a charactricstic decay length of 2.9 A. 

The polarization in two-dimensional layer is stable 
against its spontaneous reversal if 2Vyint < U, where Wint 
is the energy of interlayer interaction and U is the height 
of the potential barrier separating different orientational 
states of polarization in the layer. The multiplier two in 
the above formula relates to the worst case in which the 
considered layer is polarized antiparallel to the polariza- 
tions in neighboring layers. In our structures, the easiest 
way to reorient the polarization is to rotate it in the layer 
plane, and so the U value is equal to the energy difference 
between [110]- and [100]-polarized structures. According 
to our calculations, U — 1.84 meV for SL with n — 2 and 
1.69 meV for SL with n = 3. Comparison of U and 2Wint 
values (Table nil)) shows that the stability criterion is sat- 
isfied at n > 2. From the data presented in Tables |TT] and 
mil one can see that in superlattices with n > 3 the in- 
teraction energy between neighboring polarized layers is 
less that 10% of the ordering energy. This gives the rea- 
son to consider the ground state in these superlattices as 
an array of nearly independent ferroelectrically ordered 
planes, with quasi-two-dimensional ferroelectricity real- 
ized in each layer. 

The results obtained for (KNb03)i(KTa03)„ super- 
lattices grown on KTaOs substrates are very different 
from those obtained for stress-free superlattices. As fol- 
lows from Table IIVI for KTaOs-supported superlattices 
the ground-state structure is monoclinic for n < 4 and 
orthorhombic for n = 5-7 (data for n — 6 and 7 are 
not shown). This result can be understood as follows. 
For KNb03 grown on KTa03 substrate, the ground-state 
structure is Cm. The uniform P^j-component of polariza- 
tion docs not create macroscopic electric field in SLs, so 
that the interaction of polarizations in neighboring layers 
is weak and Px varies slowly with increasing n. On the 
contrary, the energy needed to maintain the same electric 
displacement fields normal to the layers in two materials 
increases with increasing the thickness of nonpolar layer 
n, thus weakening the instability of z-polarized unstable 
phonon. This results in fast decrease of Pz with increas- 
ing n. 

Another interesting result that follows from Table IIVI 
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TABLE IV: The ground-state structure, lattice parame- 
ters, the Berry-phase polarization (P^, Py, P^), and the 
in-plane polarization in the KNbOa layer {Px , Py) for 
(KNb03)i(KTa03)„ superlattices grown on KTaOg sub- 
strate. 



Structure Cm Cm Cm Cm Amm2 

a = 6 (A) 3.9374 3.9374 3.9374 3.9374 3.9374 

c (A) 7.9903 11.9126 15.8358 19.7591 23.6861 

Px = Py (C/m^) 0.0730 0.0605 0.0507 0.0442 0.0385 

P^ (C/m^) 0.2546 0.1769 0.1216 0.0710 0.0000 

P* = Py (C/m^) 0.1322 0.1495 0.1585 0.1646 0.1675 




is an increase of the in-plane polarization in the KNbOs 
layer with increasing n. At first sight, this result seems 
strange because the in-plane lattice parameter in SLs re- 
mains unchanged. However, if one takes into account 
the influence of Pz on the effective charge of the soft 
mode, the result can be explained. According to our cal- 
culations, in strained KNbOa grown on KTaOa substrate 
the effective charge decreased by 7.7% when the structure 
transformed from paraelectric PA/mmm to polar PAmm 
phase. Thus, when increasing n, the decrease in Pz re- 
sults in an increase of effective charge of Eu mode. This 
strengthens the dipole-dipole interaction and increases 
the instability of the soft mode. 



IV. DISCUSSION 

The critical size of ferroelectric particles and critical 
thickness of ferroelectric thin films, below which the fer- 
roelectricity vanishes, are of principle importance. First 
studies of ferroelectric nanoparticles indicated that the 
ferroelectricity vanishes when the particle diameter is 
about 100 A. However, the experiments on ultrathin fer- 
roelectric filmsSE and theoretical calculations from first 
principles^"— have shown that a non-zero polarization 
normal to the film surface is retained in films with a thick- 
ness of 3-6 unit cells if the depolarizing field is compen- 
sated. In the PVDF-TrFE copolymer films the switch- 
able polarization was observed in films with a thickness 
of 10 A (two monolayers).— The ferroelectric state with 
the polarization parallel to the layers was predicted for 
BaTiOa and PbTiOs films with a thickness of 3 unit 
cells^i^ and PZT films with a thickness of one unit 
celli^ Unstable ferroelectric mode with the polarization 
in the film layer was predicted for films of PhB'^^^By^Oz 
solid solutions with a thickness of one unit cell.— In all 
these cases two-dimensional ferroelectricity was a result 
of decreasing the physical size of the samples in one di- 
rection. In contrast, in this work, quasi-two-dimensional 
ferroelectricity appears as a result of specific interactions 
in superlattices and can be obtained in structures of ar- 
bitrary thickness. 



FIG. 3: Phonon dispersion curves for paraelectric PA/mmm 
phase of (KNbOs)! (KTa03)3 superlattice. 



The polarization profiles obtained in this work agree 
with results of atomistic modeling of KNbOa /KTaOa su- 
perlatticesJ^i^ According to the data of these papers, 
the polarization parallel to the layers decreased by 3-4 
times when shifting by one unit cell from the hetero- 
junction (in contrast to the behavior of the polarization 
normal to the layers) . Our conclusion about the instabil- 
ity of phonon in the paraelectric phase of SLs with 
n < 3 agrees with results of Ref. [l^ only partially: in the 
latter paper stable polar phase with P4mm symmetry 
was obtained only in SLs with n = 1 and 2. The insta- 
bility of the Eu mode in SL with n — 7 agrees with that 
observed in Ref. [2l|, but in this work it is much stronger. 

However, in spite of qualitative agreement of results of 
this and previous works, the physical conclusions that 
we make are very different. First, in all stress- free 
KNbOs/KTaOs superlattices the PAmm phase is not the 
ground-state structure, as well as in KTaOs-supported 
superlattices with n > 5. The reason for this is a ten- 
dency of polarization in SLs to incline towards the layer 
plane. This tendency was revealed in a number of fer- 
roelectric superlattices^'^ and was shown to be a way to 
reduce the electrostatic and mechanical energies in them. 

Second, the existence of stable polarization in KNbOa 
layers with a thickness of one unit cell and very inho- 
mogeneous distribution of polarization, which is accom- 
panied by exponential decay of the interaction energy 
between neighboring layers with increasing the inter- 
layer distance, result in a new, previously unknown fea- 
ture of the ground-state structure of (KNb03)i(KTa03)„ 
superlattice — the formation of an array of nearly inde- 
pendent ferroelectrically ordered planes in the bulk of 
the superlattice. 

Phonon dispersion curves calculated using the interpo- 
lation technique from phonon frequencies at F, X , M , Z, 
R, and A points of the Brillouin zone for PA/mmm phase 
of (KNb03)i(KTa03)„ superlattices enables to clarify 
the origin of unstable phonon modes at the boundary of 
the Brillouin zone in these SLsj-^ " The dispersion curves 
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for these modes are shown in Fig.[3l Analysis of the eigen- 
vectors shows that the doubly degenerate unstable mode 
along T-Z direction is polarized along [100] and [010] di- 
rections while nondegenerate unstable modes along Z~R, 
R~X, and X~T directions are polarized along [010] di- 
rection. The character of atomic displacements for these 
modes is fully consistent with the ferroelectric instability 
of ...-0-Nb-O-... chains first observed in KNbOsf^a but 
in SLs the chains are formed only in [100] and [010] di- 
rections. A nearly zero dispersion of the unstable mode 
along T-Z direction means that the interlayer interaction 
is weak, and the instability against the formation of fcr- 
roelectrically ordered planes does not depend on whether 
neighboring planes are ordered ferroelectrically or anti- 
ferroelectrically. This supports the idea about indepen- 
dence of ferroelectrically ordered planes in SLs with thick 
KTaOa layers. 

The obtained results demonstrate that in complex sys- 
tems in some cases the polarization calculated using the 
Berry phase method cannot be interpreted as a bulk 
property. The variations in polarization in the unit cell 
can be so large that the ferroelectric system can acquire 
the properties of two-dimensional system. To avoid inter- 
action between polarized layers, the polarization should 
be parallel to the layers (a non-zero normal component 
results in appearance of electric fields in the structure, 
and 2D behavior transforms to 3D one). This means 
that to create arrays of nearly independent ferroelectri- 
cally ordered planes in [001]-oriented superlattices, we 
need SLs with Pmm2 or Amm2 ground-state structure. 

Complex ferroelectric/antiferroelectric behavior has 
been observed recently in Ruddlesden-Popper (RP) su- 
perlattices,—!^ in which appearance of Goldstone-like ro- 
tational mode of the in-plane polarization is possible4^ In 
stress- free BaTiOa/SrO superlattices, the in-plane polar- 
ization appeared in SLs with the BaTiOs slab thickness of 
n > 4 unit cells; when these SLs were epitaxially grown 
on SrTiOa substrates, the polar structure transformed 
to antiferroelectric (domain-like) structure in which out- 
of-plane polarization alternated along [100] or [110] di- 
rections in each of BaTiOs slabs. The antiferroelectric 
structure was studied in more detail in BaTiOs/BaO su- 
perlattices grown on SrTiOa substrate. "^^ The instability 
at X point of the Brillouin zone appeared in structures 
at a thickness of BaTiOa slabs of 2-3 unit cells while 
the zone-center ferroelectric instability appeared only at 
n = 4-5 unit cells. It was supposed that a huge depo- 
larizing field appears in SLs because of the difficulty to 
polarize BaO layer, and this strongly harden the fer- 
roelectric mode at the F point. At the same time, the X 
phonon mode arising from the same ferroelectric insta- 
bility of ...-0-Ti-O-... chains does not produce macro- 
scopic depolarizing field, and so the influence of the BaO 
layer on the phonon frequency is weaker. This explains 
why in RP superlattices the antiferroelectric phase can 
appear before the ferroelectric one. 

The superlattices studied in this work differ markedly 
from RP ones^i*^ because in our case the ferroelectric 



layers are separated by the layers of incipient ferroelec- 
tric, which has much higher dielectric constant compared 
to SrO (BaO). In contrast to the RP superlattices, our 
SLs do not exhibit antiferrodistortive instability. In addi- 
tion, the thickness of the KNbOa layer which is necessary 
to induce ferroelectricity in our SLs is only one unit cell. 
To test the influence of the thickness of KNbOa layer on 
the appearance of aforementioned antiferroelectric struc- 
ture, we performed additional calculations for stress-free 
SLs in which KNbOa layers with a thickness of 3 and 
5 unit cells were separated with a thick (3 unit cells) 
KTaOa layer. 

The calculations showed that all three 
(KNb03)„(KTa03)3 SLs with n = 1, 3, and 5 studied in 
this work are unstable against the appearance of antifer- 
roelectric structure with out-of-plane displacements of 
ions described by X phonon [Pmrna space group), but 
in all cases the energies of corresponding equilibrium 
structures are higher than the energies of the ferroelec- 
tric Amm2 phases. This means that antiferroelectric 
Pmma structure is not the ground state of stress-free 
KNbOs/KTaOs superlattice. Moreover, the calculation 
of phonon frequencies at the F point for Pmma phase of 
(KNb03)i(KTa03)i SL showed that this antiferroelec- 
tric phase exhibits the ferroelectric instability against 
the appearance of in-plane polarization. Therefore one 
can conclude that like for RP superlattices, the antifer- 
roelectric phase of KNbOa/KTaOa superlattices with 
out-of-plane displacements of ions is indeed more stable 
compared to the ferroelectric phase with out-of-plane 
polarization, but nevertheless the ferroelectric Amm2 
phase with non-zero in-plane polarization has the lowest 
energy. In RP superlattices studied so far, the in-plane 
polarization was suppressed by epitaxial strain from the 
substrate, and the antiferroelectric phase was the most 
stable one. 

We discuss now possible applications of the ar- 
rays of nearly independent ferroelectrically ordered 
planes. At first sight, the fact that the polarization 
in KNbOa/KTaOa superlattices is oriented in the layer 
plane can be regarded as not intriguing. So far, when 
discussing the ferroelectricity in ultrathin films, main at- 
tention was paid to films in which the polarization is nor- 
mal to the film surface as these films are more suitable 
for different applications. However, the possibility of for- 
mation of an array of nearly independent ferroelectrically 
ordered planes opens new opportunities in design of elec- 
tronic devices. These structures can be used as a medium 
for three-dimensional information recording. Accepting 
that the lateral size of a ferroelectric domain that has a 
long-time stability against spontaneous reversal of polar- 
ization is 250 A,^* for the interlayer distance of 16 A (the 
period of SL with n — 3) the estimated 3D information 
storage density can be as high as ~10^® bit/cm'^. This 
value exceeds the density achieved in modern optical 3D 
recording media by six orders of magnitude. 

One can add that the two-component order parameter, 
which admits the rotation of polarization by 90° in each 
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of independent KNbOa layers, offers one more opportu- 
nity. In ferroelectric RAM devices with polarization nor- 
mal to the film surface the leakage currents and conduc- 
tivity of films are main effects preventing the realization 
of non-destructive read-out. The systems admitting the 
rotation of polarization by 90° have a feature that en- 
ables to overcome this disadvantage — a non-destructive 
read-out of polarization in them can be realized using the 
anisotropy of the dielectric constant in the layer plane. 
Our calculations predict that in (KNb03)i(KTa03)„ su- 
perlattice with n — 2 the dielectric constants parallel 
and perpendicular to the polarization are 142 and 188, 
whereas in the SL with n — 3 they are 225 and 281. 

V. CONCLUSIONS 

In this work, the density-functional theory was used to 
calculate the properties of (KNb03)i(KTa03)„ superlat- 
tices with n = 1-7. It was shown that the ferroelectric or- 



dering is retained in KNb03 layers of minimal thickness 
(one unit cell) whereas the interaction energy between 
neighboring polarized layers decreases exponentially with 
increasing parameter n. Array of nearly independent fer- 
roelectrically ordered planes with the polarization lying 
in the layer plane was shown to be the ground-state struc- 
ture for SLs with n > 2. Using of these structures as a 
medium for three-dimensional (3D) information record- 
ing enables to attain the information storage density of 
IQis bit/cm3. 

The calculations presented in this work were performed 
on the laboratory computer cluster (16 cores) and SKIF- 
MGU "Chebyshev" supercomputer at Moscow State Uni- 
versity. 
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